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Abstract 

Background: Cliagas disease, caused by Trypanosoma cruzi, is a geograpliically widespread antliropozoonosis tliat 
is considered a major public health problem in Latin America. Because this parasite presents high genetic variability, 
a nomenclature has been adopted to classify the parasite into six discrete typing units (DTUs): Tel, Tell, Tclll, TcIV, 
TcV, and TcVI, which present different eco-epidemiological, clinical, and geographic associations. Currently, the 
available genotyping methods present a series of drawbacks that implies the need for developing new methods for 
characterizing T. cruzi DTU's. The aim of this work was to genotype reference populations from T. cruzi by means of 
a High-Resolution Melting (HRM) genotyping assay. 

To genotype the DTUs of 38 strains and 14 reference clones of T. cruzi from diverse sources, real-time PCR (qPCR) 
was coupled to high-resolution melting (HRM) based on the amplification of two molecular markers — the divergent 
domain of the 24 sa rRNA gene and the intergenic region of the mini-exon gene. 

Findings: Amplification of the mini-exon gene allowed the genotyping of three distinct groups: Tel, Tell- TcIV- TcV, 
and Tclll-TcVI, while amplification of the 24sa gene generated non-overlapping melting temperature ranges for 
each DTU that were used to identify the groups in the six existing DTUs of Trypanosoma cruzi. 

Conclusions: The proposed genotyping assay allowed discrimination of the six genetic groups by obtaining 
specific melting curves for each DTU. The application of this technique is proposed because of its specificity, 
sensitivity, high performance, and low cost compared with other previously described characterization methods. 
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Findings 

Chagas disease is an anthropozoonosis caused by the 
parasite Trypanosoma cruzi and affects an estimated 10 
milUon people, posing a serious public health problem 
[1]. Trypanosoma cruzi has a high genetic variability with 
evidence of six Discrete Typing Units (DTUs): T, cruzi I 
(Tcl), T cruzi II (TcII), T cruzi III (TcIII), T cruzi IV 
(TcIV), T cruzi V (TcV), and T cruzi VI (TcVI), each of 
which presents characteristics based on the geographic 
distribution, eco-epidemiological associations, and clin- 
ical manifestations of the disease [2]. Currently, several 
methods have been reported for characterizing and 
genotyping Trypanosoma cruzi populations, including 
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the PCR amplification of the SL-IR (spliced-leader 
intergenic region of the mini-exon gene), 24sa, and 18s 
rDNA gene regions. The comparison criterion for this 
technique is based on the size differences in the resulting 
amplicons [3,4]. However, the characterization obtained 
using this technique is limited by the absence of ampli- 
fication products or the need for multiple amplification 
procedures [2,5]. Another proposed methodology is 
based on genetic polymorphism analysis using multilocus 
PCR-RFLP; however, this technique requires the develop- 
ment of combinations and multiple assays, increasing its 
complexity and limiting its use [6,7]. In 2009, Llewellyn and 
colleagues implemented an approach based on the 
Multilocus Microsatellite Typing (MLMT) of a panel of 48 
polymorphic microsatellite markers. This technique is a 
quick and feasible alternative but requires the inclusion of 
sequencing procedures that lead to increased costs [8]. A 
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Table 1 Biological and geographical origin of the 36 
reference strains and 18 reference clones of Trypanosoma 
cruzi previously genotyped by MLEE, MLST and MLMT 
employed in the HRM analyses 



Table 1 Biological and geographical origin of the 36 
reference strains and 18 reference clones of Trypanosoma 
cruzi previously genotyped by MLEE, MLST and MLMT 
employed in the HRM analyses (Continued) 



Tel 



Tell 



DTU Strain/clone 

XI 0610 
92101601 
Cutia ell* 
DM28 
Jr CI4* 
xlOl 
SaxplS 
Chilee22* 
lOr 26 
c8 
ERA 

SP104eir 
P209el93* 
0PS21elll* 
MASeir 
IWel4* 
Pot 7b 
CBB 
Pot 7a 
Chaeo 23 
EBel2r 
Tul8 
ESM 
85/847 
M6421 
SLDNl 
Ua2 
SABPI9 
M5631el5* 
StelOr 
92122 
TelV CANIII 
TV 

DOG THEIS 
92.80 
Chaeo 9 
Para6 

PPAH179 CL5* 
Chaeo 2 CB 
AACF2 
Para4 
se43 



Telll 



TeV 



Geographical origin 

Guarieo, Venezuela 
Georgia, USA 
Espiritu Santo, Brazil 
Carabobo, Venezuela 
Anzoategui, Venezuela 
Para, Brazil 
Maje, Peru 
FlordeValle, Chile 
Santa Cruz, Bolivia 
La Paz, Bolivia 
Anzoategui, Venezuela 
Region IV, Chile 
Suere, Bolivia 
Cojedes, Venezuela 
Minas Gerais, Brazil 
Cuneumen, Chile 
San Martin, Paraguay 
Tulahuen, Chile 
San Martin, Paraguay 
Chaeo, Paraguay 
Boyaea, Colombia 
Tupiza, Bolivia 
Sao Felipe, Brazil 
Alto Beni, Bolivia 
Belem, Brazil 
Casanare, Colombia 
Suere, Bolivia 
Vitor, Peru 
Marajo, Brazil 
Georgia, USA 
Georgia, USA 
Belem, Brazil 
Boyaea, Colombia 
Oklahoma, USA 
Santa Cruz, Bolivia 
Chaeo, Paraguay 
Chaeo, Paraguay 
Chaeo, Argentina 
Chaeo, Paraguay 
Casanare, Colombia 
Paraguari, Paraguay 
Santa Cruz, Bolivia 



Biological origin 

Homo sapiens 
Didelphis marsupialis 
Dosyprocto oguti 
Didelphis marsupialis 
Homo sapiens 
Homo sapiens 
Homo sapiens 
Triatoma spinolai 
Aotus sp. 

Triatoma infestans 
Homo sapiens 
Mepraia spinolai 
Homo sapiens 
Homo sapiens 
Homo sapiens 
Homo sapiens 
Triatoma infestans 
Homo sapiens 
Triatoma infestans 
Triatoma infestans 
Homo sapiens 
Triatoma infestans 
Homo sapiens 
Dasypus novemcinctus 
Homo sapiens 
Dasypus novemcintus 
Dasypus novemcintus 
Triatoma infestans 
Dasypus novemcinctus 
Procyon lotor 
Procyon lotor 
Homo sapiens 
Triatoma venosa 
Canis familiaris 
Homo sapiens 
Triatoma infestans 
Triatoma infestans 
Homo sapiens 
Triatoma infestans 
Canis familiaris 
Triatoma infestans 
Triatoma infestans 



TeVI 



VinehlOl 


Limari, Chile 


Triatoma infestans 


MN el2* 


Region IV, Chile 


Homo sapiens 


BUG2148ell* 


Rio Grande do Sul, Brazil 


Triatoma infestans 


SC43eir 


Santa Cruz, Bolivia 


Triatoma infestans 


S03el5* 


Potosi, Bolivia 


Triatoma infestans 


LhVa 


Chaeo, Argentina 


Triatoma infestans 


VfrA 


Franeia, Chile 


Triatoma infestans 


Tula CI2* 


Tulahuen, Chile 


Homo sapiens 


Rp540 


Casanare, Colombia 


Rhodnius prolixus 


p251 


Coehabamba, Bolivia 


Homo sapiens 


CLBRENER 


Rio Grande do Sul, Brazil 


Triatoma infestans 


P63ell* 


Makthlawaiya, Paraguay 


Triatoma infestans 



^Biological Clones. 

similar initiative was proposed to develop a Multilocus 
Sequence Typing (MLST) scheme based on the analysis of 
a panel of four molecular targets: Met-III, RB19, TcGPXII, 
and DHFR-TS; however, the development of this process 
requires both bioinformatic analysis and sequencing steps 
[9]. Recently, the analysis of the polymorphisms of the 24sa 
rDNA gene, the subunit II of the cytochrome oxidase, and 
four microsatellite loci has been shown to be an alternative 
method [10]. This technique is highly efficient but requires 
multiple amplification steps that can lead to an increased 
risk of contamination. 

Thus, a low-cost genotyping technique with reliable 
and robust results for the identification of the six genetic 
groups in a single assay is needed to facilitate the 
characterization of the circulating strains in a specific 
geographic area. A technique such as High-Resolution 
Melting (HRM) offers a viable alternative for genotyping 
T. cruzi given its robustness, speed, cost, and the elimin- 
ation of sequencing or hybridization procedures [11]. This 
molecular analysis has recently been used to genotype 
different populations of parasitic microorganisms [12-15]. 
The aim of this study was to genotype 36 reference strains 
and 18 reference clones of Trypanosoma cruzi previously 
genotyped by different molecular markers from different 
biological and geographical origins via HRM. 

We employed 36 strains and 18 reference clones of 
Trypanosoma cruzi previously characterized via MLST, 
MLEE, MLMT, and RFLP-PCR (Table 1). Additionally, 
in order to test specificity we included DNA from reference 
strains of T, rangeli and Leishmania panamensis, A 
real-time PGR assay coupled with HRM analysis was 
developed to target SL-IR and the 24sa rRNA gene using 
a 7500 Fast Real-Time PGR System (Applied Biosystems, 
Inc., GA, USA). These genes were selected because they 
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TcVl 
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Tel 



Tell 

0,23 



Telll 

0,16 



TcIV 

0,08 



TeV 

0,23 



TcVI 
0,41 



0,26 

Discrete Typing Units (DTU's) ^ 

Figure 1 Melting Temperature profiles and scatter plots of the HRM genotyping assay. A. Representative Normalized melting curves of the 
Trypanosoma cruzi Discrete Typing Units B. Scatter-plot showing overlapping and non-overlapping changes in the Tm values of the curves 
obtained by each DTU; the values above the boxes indicate the number of strains/clones employed per DTU. 



had previously been used as molecular markers for the 
PGR genotyping of T. cruzi [4,16]. The amplification 
reactions were carried out in a total volume of 21 (iL. The 
reaction mixture consisted of Ix MeltDoctor HRM Master 
Mix (Applied Biosystems, Inc., CA, USA); a 5 (iM solution 
of each of the primers TCC {S'CCC CCC TCC CAG GCC 
ACA CTG 3) and TC2 (5'CCT GCA GGC ACA cGT GTG 
TGT G 3) for the amplification of SL-IR and 24sa, D71 
(5'AAG GTG CGT CGA CAG TGT GG 3) and D72 
(5'TTT TCA GAA TGG CCG AAC ACT 3); 6.6 (iL of 
water; and 20 ng/(iL DNA template. The cycle condi- 
tions for the qPCR were as follows: 95°C for 10 minutes 
(1 cycle), followed by amplification for 40 cycles at 95°C 



for 15 seconds (denaturation) and 60°C for 1 minute. 
The amplicon dissociation was immediately initiated by a 
melting step. The thermal profile consisted of denaturation 
at 95°C for 10 seconds, 60°C for 1 minute (annealing), 95°C 
for 30 seconds (high-resolution melting), and a final 
annealing stage at 60°C for 15 seconds. The melting curve 
analysis was carried out by supplying raw fluorescence (F) 
vs. temperature (T) values or already transformed -dF/dT 
values. The raw value,differentiation was performed and 
the first derivative curve is displayed with the melting 
peaks. Peaks and Tm values were automatically identified, 
peak areas calculated and peaks flagged as negative' if the 
peak area is below the user-defined threshold value. During 
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this process, the PGR amplicons were denatured prior 
to developing the melting curves, and the changes in 
fluorescence with respect to the changes in temperature 
(dF/dT) were recorded by the instrument with a ramp 
of 0.5°C/seg. Each DNA sample was analyzed in duplicate. 
Finally, the normalized and derivative melting curve 
profiles were obtained. The analysis was performed 
using a blind methodology to reduce the possibiUty of 
bias. To evaluate the reproducibility of the assay, the 
genotyping technique was performed on five independent 
days under the same conditions. The specificity was 
assessed by evaluating the concordance between the 
pre-established genotypes of the DNA samples used and 
the genotypes characterized by HRM. The consistency of 
the method was determined by correlating the previously 
quantified reproducibility and specificity [17]. Finally, a 
randomization and subsequent non-parametric Kruskal 
Wallis test (p < 0.05) were conducted to detect any signifi- 
cant differences between the melting temperature data 
for each DTU. The analyses were performed using the 
R platform and GraphPad software package. 

A specific melting temperature was obtained for each 
of the reference strains and clones included in this study, 
both for amplifying SL-IR and for the 24sa. The plots 
were created for the melting temperature distribution of 
each DTU (Figure 1). The 24sa gene exhibited higher 
discriminatory power than the SL-IR gene, producing 
non-overlapping ranges of the melting temperature for 
each DTU. The normalized fluorescence curve and the 
derivative melting curve yielded a different and easily 
distinguishable profile for each DTU. The representative 
profiles were collected for the derivative, the normalized 
melting curves for the SL-IR gene and the 24sa amplicons 
(Figure 1). Finally, an algorithm was constructed (Figure 2) 
which allowed the discrimination of the DTUs after the 
amplification of the SL-IR gene and the subsequent ampli- 
fication of the 24sa gene. To validate the method, three 



parameters were determined based on the methodology of 
Moncayo and Luquetti (1990) [17]; the reproducibility 
was determined for both the SL-IR gene amplification and 
the 24sa rRNA and yielded a value of 100% for both tests. 
The specificity and consistency values were obtained only 
for the amplification of the 24sa rRNA, which was the 
only marker capable of generating non-overlapping melting 
temperature ranges that was specific for the six DTUs and 
the DNAs from T. rangeli and L panamensis. 

The use of the HRM methodology allowed the adequate 
and congruent discrimination of the T, cruzi genotypes for 
the strains and reference clones used. Many problems, 
such as the need for sequencing, the degree of complexity, 
the design of multiple pairs of primers, and genotyping by 
the absence of a band, among others, have been reported 
with the existing genotyping methods. Thus, new alterna- 
tives are needed to simplif)^ the process and minimize the 
difficulties. For example, methods based on the absence 
of an amplification product are considered problematic 
because no actual evidence is generated that may relate 
to the presence of a specific DTU [2]. To obtain a correct 
identification, other genotyping methods for which high 
efficiency, sensitivity, and robustness have been reported 
require multiple amplifications of a single sample using 
different markers or coupled with RAPD s and RFLP 
analysis, which can significantly increase the economic cost 
and lead to sample contamination [18]. The HRM tech- 
nique does not present any of the aforementioned difficul- 
ties and, given the as-obtained results, could be considered 
a potential alternative strategy for genotyping T, cruzL 

The markers were selected based on two criteria: the 
fact that both the SL-IR gene and the 24sa rRNA have 
previously been used to characterize T, cruzi and because 
they exhibit a high number of copies within the parasite 
genome, which increases the probability of detection 
[6,19,20]. The results of this study indicated that the SL-IR 
gene did not aUow the absolute discrimination of the 
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SL-IR 



(TCC/ TC2) 



24sa 



(D71/D72) 



Tc I 



Tc I 



TmValue 



86.010.4 85.5±0.1 85.4±0.4 84.3±0.3 84.6±0.2 84.1 ±0.3 

Z3 ljD l □ [ 



Tclll/TcVI 



Tell / TcIV /TcV 



80.4±0.3 76.310.3 79.210.2 79.810.2 78.710.2 77.510.4 



Tc III 



Tc VI 



Tc II 



Tc IV 



TcV 



Figure 2 Feasible and accurate algoritiim siiowing tiie use of HRM assay for thie discrimination of Trypanosoma cruzi DTU's. 
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six DTUs in a single assay, a fact that may be due to the 
recently confirmed variation in the copy number of the 
gene and the chromosomal rearrangement at the level of 
the different strains [21]. To examine the reproducibility 
of the method, the assays were conducted in duplicate, 
and the tests were repeated over several days, maintaining 
the same test conditions, reagent quantities, and DNA 
concentrations. The technique generated a reproducibility 
value of 100% for both molecular markers. This result is 
of great significance because genotyping techniques must 
be reproducible and robust to avoid producing conflicting 
results and to accurately establish the genotype of each 
sample analyzed [22] . 

Conclusion 

The use of a real-time PGR assay coupled with HRM 
analysis is proposed for genotyping the DTUs of the 
T. cruzi population for the development of research 
projects that require genotyping of the circulating 
strains because it provides an easy, quick, and low 
(1.5 USD per assay compared to 2.5 USD for genotyping) 
alternative compared to the existing methods. Among the 
advantages over other currently used methods are its speed 
and simplicity (given that this method does not require gel 
electrophoresis), the reduction of the processing time, and 
the specificity and sensitivity of the method [14]. Therefore, 
HRM coupled to qPCR provides a cost-effective alternative 
for the future genotyping and identification of the six 
DTU s of T. cruzi which fulfills one of the main research 
targets in Chagas disease. 
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